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In this work, we illustrate how a Jordan- Wigner transformation combined with symmetry con- 
siderations enables a direct solution of Kitaev's model on the honeycomb lattice. We (i) express 
the p-wave type fermionic ground states of this system in terms of the original spins, (ii) adduce 
that symmetry alone dictates the existence of string and planar brane type correlators and their 
composites, (iii) compute the value of such non-local correlators by employing the Jordan- Wigner 
transformation, (iv) affirm that the spectrum is inconsequential to the existence of topological quan- 
tum order and that such information is encoded in the states themselves, and (v) express the the 
local symmetries of Kitaev's model and the anyonic character of the excitations in terms of fermions. 

PACS numbers: 



I. INTRODUCTION 

Topological quantum order (TQO) [l[ Q is a new 
paradigm that lies beyond the realm of Landau's theory 
[3| [J] of local order parameters. TQO is intuitively asso- 
ciated with insensitivity to local perturbations: the order 
is topological. As such, TQO cannot be described by local 
order parameters. A quintessential example of a system 
with TOO is Kitaev's model on the honeycomb lattice 
@1 [1] 0- br this article, we will study this model. Our 
new central results are (a) an explicit form for the ground 
states in real space which extends and complements the 
works of @ and @ as well as (b) our finding of non-local 
correlations [two dimensional string or brane type correla- 
tors] in this system. We further show how the (c) known 
anyonic excitations of Kitaev's model can be examined 
anew by using a direct Jordan- Wigner transformation. 
Beyond providing a direct solution which highlights cer- 
tain previously overlooked aspects, our results will flesh 
out some of the more general ideas @ , [1] regarding the 
general character of TQO. 



II. OUTLINE 

This article is organized as follows: In section IIII1 
we review the fermionization of the Kitaev model on a 
hexagonal lattice 0, Qlj. The original model of spins on 
a hexagonal lattice is mapped to a model of p-wave BCS 
model with sitc-dcpcndcnt chemical potential for spinlcss 
fermions on a square lattice. In section lPvT, we discuss the 
fermionic representation of symmetries embodied in the 
model. It is shown that the local conserved quantities 
on the plaquettes studied by Kitaev § are equivalent to 
the conserved bond quantities in the fermionic represen- 
tation, up to a gauge fixing. In section [Vj the ground 
state configuration, which is vortex free in spin represen- 
tation and with a uniform chemical potential in fermionic 
representation, is exactly solved in fermionic representa- 



tion. The spin basis form of the ground state is also 
studied. It is shown that the ground state can be writ- 
ten as a projection of a given reference state over a sec- 
tor of Hilbert space that has a chosen set of topological 
numbers which are similar to, yet slightly more compli- 
cated than, other systems which exhibit TQO such as the 
Rokhsar-Kivelson dimer model[ll[ and the Kitaev toric 
code model@. The phase diagram and Bogliubov exci- 
tations are also obtained. In Section (|VI[) , we show how 
the known results about the anyonic vortex excitation 
in gapped state can be easily studied in the fermionic 
representation. In section IVII1 a symmetry argument is 
put forward concerning the vanishing correlation func- 
tions, which is also recently derived anew by a Majo- 
rana fermionization construction [l2j. We also show that 
string correlators naturally appear in Kitaev's model. A 
brief summary in scction lVIIII concludcs this work. In the 
appendices, we provide technical details concerning the 
determination of the ground state and review Elitzur's 
theorem. 



III. FERMIONIZATION 

We begin with a Fermionization of the Kitaev model 
on the hexagonal lattice Q [|| which is defined by the 
following S = 1/2 Hamiltonian 

H = -J x J2 °R°R'- J y E 

x — bonds y — bonds 

z — bonds 

with R and R' lattice sites. This system can be fermion- 
ized by a Jordan- Wigner transformation]^, [l(| ■ This one- 
dimensional fermionization is made vivid by deforming 
the hexagonal lattice into a "brick-wall lattice" which is 
topologically equivalent to it on which we may perform 
a one dimensional Jordan- Wigner transformation. The 
schematics are shown in FigJT] In the up and coming, we 
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will mark all sites by "white" or "black" (w/b)in order 
to denote to which sublattice which they belong to. The 
distance between two nearest-neighboring sites on this 
lattice will be set to unity. 
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FIG. 2: Schematics of the contour for the Jordan- Wigner 
transformation that we employ in the deformed hexagonal 
lattice- see text and Eq.© in particular. 



FIG. 1: Deformed hexagonal lattice and three types of bonds. 

Throughout, we will consider the system with open 
boundary conditions unless stated otherwise. The di- 
agonal directions e x and e y shown in Fig.JT]) will be of 
paramount importance in our final solution. Let us de- 
note by the Cartesian coordinates of each site R on 
the lattice of Fig. [TJ Let us next consider the Jordan- 
Wigner transformation defined by a simple one dimen- 
sional contour which threads the entire lattice [See Fig[2]: 



r£ = 2 



nn< 



i'j' 



j'<j i> 



n 



2c!„-c, 



1. 



This path goes through each lattice site exactly once as 
shown in Fig.©. In Eq.©, cr+ = (a x + ia y ) is twice the 
spin raising operator at a given site- hence the factor of 
two. 

The Kitaev model of Eq. |T]) now becomes 

h = Jx y, ( ct - c U ct + c ) 6 

x — bonds 

•y E 

y — bonds 



Jy 

-J, 



(ct + c) 6 (ct 



i) b ( 2c t c 



(3) 



z — bonds 



Henceforth, the subscripts b and w will denote the white 
and black sites of a bond as illustrated in Fig[T] 
Let us next introduce the Majorana fermions 

(c - S) w ji B n 



A w = (c-c 1 ) w /i B w = (c+c r ), 
for the white sites and 

B b = (c-J) b /i A b = (c + J) b 



(4) 
(5) 



for the black sites. With all of these transformations in 
tow, the Hamiltonian now reads 



H = -i 



E JxAw-A-b Y^ JyAt,A v 

y — bonds 



x — bonds 



-J, Jz{BA) b {BA) l 



(6) 



-bonds 



It is easy to see that BB along the z-bond is a conserved 
quantity [Io| . Thus, the Z 2 operator 



a r — iB b B w , 



(7) 



with r the coordinate of the midpoint of the bond con- 
necting the black and red sites, is fixed for each vertical 
bond. The Hamiltonian of Eq.|T|) now reads 



H{{a}) = -i 



y — bonds 



x — bonds 



iJ z OL r AbA v 



(8) 



z — bonds 



Here, r denote the centers of the vertical bonds. In 
Section {TV}, We will show that {a r } are intimately re- 
lated to the local symmetries present in Kitacv's model 
of Eq. |T]) . This identification, combined with Reflection 
Positivity arguments, Q will allow us to infer that, up to 
[d = 1 [l3|) symmetry operations, a r = 1 for all r. The 
ground state does not contain any "vortices" which are 
marked by one dimensional in the Ising variables {a r } 
along a row. Similar Reflection Positivity arguments re- 
garding the absence of vortices in other systems and a 
bound on the energy penalties that they entail are. e.g., 
given in [hH ]. This, in turn, will allow us to explicitly 
diagonalize the Hamiltonian. 
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IV. FERMIONIC REPRESENTATION OF 
LOCAL SYMMETRIES 

As shown by KitaevjBj, the Hamiltonian of Eq.fl]) has 
one conserved quantity for each plaquctte (or hexagon) 
h, 



4 e 4 a 4 e 4 e 4 e 4 



t y z x y z x 

1 h = (J \ w (7 2b (T ?,w (J Ab (T bw (J Qb- 



(9) 



Here, 1 — 6 denote the sites of a given plaquette, as illus- 
trated in the plaquette of the lower corner in Figfl] In 
the subscripts, we also label the (w/b) sublattices of the 
sites. The conserved quantity Ih on plaquette is equiva- 
lent to the bond conserved quantity a defined in previous 
section, up to a gauge fixing as we shall show below. This 
can be shown by fermionizing Ih using the transformation 
we introduced in previous section. Let us first fcrmionizc 
the product of first three spins 



y z x 
a \w a 2b a -iw 



l lw (J 2b CT lw (J 2b 

= iB lw B 3w . 



,t 



' 3 <r 



3w 



Similarly, 



Therefore, 



^tb a lm^lb = iB Ab B t 



<)/>■ 



a 3 40!i 6 . 



(10) 



(11) 



(12) 



This model has an extra degeneracy that links differ- 
ence sectors parameterized by different sets of a's. The 
vortex [or anyon] variables are the product of two consec- 
utive Ising bond variables. In other words, anyons (Ih) 
are none other than domain walls in the Ising variables 
({a r }) that our system contains. Inverting all of the val- 
ues of a r for all sites r which lie along a horizontal row 
leaves the system unchanged. Physically, effecting the 
transformation a r — > — a r for all bonds r along a row 
does not change the vorticity content of the system: all 
domain walls along the chain remain invariant (and as we 
show so does the spectrum). In the notation of [|[ and 
[T3 |. this corresponds to a d = 1 dimensional operator 
[as it involves spin operators on a (d = 1 dimensional) 
line]. For simplicity, let us divide a into different subsets: 
{a} = L>i[a]i. Here, [a]i denotes the set of bonds that 
are connected to white sites of the i-the horizontal line 
of the brick-wall lattice. We next explicitly write down 
these d = 1 symmetry operators (see [8| , [13| for a defini- 
tion of d dimensional symmetry operators) in both their 
fermionic and original spin language form. Towards this 
end, we construct the unitary operator U™ 



(13) 



In Eq. (fT5|) . n is a site index. The product is taken over 
all sites in the rows j > i. We notice that Uf effectively 



f i = l 



+■ i = l 



f i = l 



+■ 1 = 1 +■ i = l + i = l + i = l + 1 = 1 + ] 



f i = l 



■i- i = l +■ i = l f 1 = 1 f 1 = 1 f 1 = 1 
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(b) 

FIG. 3: (a) Vortex free configuration. (b)Reversing / in a 
plaquette corresponds to flip the chemical potential along a 
string. 



reverses the sign of [a]i while leaving others untouched 



= H(..., [a],_i, -[a]i, [a] i+ i, ...). 



(14) 



Therefore, we have extra freedom to fix one a in each sub- 
set [a]i. This degree of freedom is closely related to the 
string-like conserved quantity in the original spin model, 
namely, 



P, 



n 



(15) 



which rotates all spins on j-th line by 180 degree around 
z. Under the gauge fix in which one a is fixed to +1 in 
each [a]i, the set of Kitaev conserved quantities {Ih} is 
equivalent to {a}. 



V. GROUND STATE 
A. Diagonalization 

Armed with a physical meaning of the bond variables 
a from Section( |IV[) . we now proceed to solve the problem 
posed by Eq.® and diagonalize the problem. 

It is interesting to introduce a fermion on the z-bond 



{A w +iA b )/2, S = {A w -iA b )/2, 



(16) 



where A w and A b arc the Majorana fcrmions on the white 
and black site of a given z-bond. We thus have a model 
for fermions on a square lattice with a site dependent 
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chemical potential 



r 

+ Jy j2(dt+d r ) (4+6,-* 

r 



r+e. 



(17) 



J, 



J, 



The unit vector e y connects two z-bonds and crosses a y- 
bond,see Fig. ([I]). A similar definition holds for e x . Note 
that 



[a r ,d r ] = [otr,4] = 0- 



(18) 



For sufficiently large systems, the ground state config- 
urations are bulk vortex-free configurations [f| in which 
Ih = 1 for all plaquettes (hexagons) h. The ground states 
of the fcrmionic problem of Eq. I|17p has a r = 1 everwhere 
[and all sectors of {a r } related to it by the d = 1 oper- 
ation of Eqs. (|13ll5p ] which corresponds to //, = 1 for 
all plaquettes h. As shown in Fig. ([3]), reversing Ih to 
Ih = — 1 leads to an inversion of the chemical potentials 
{a r } along a horizontal string. 

The exact solution for ground state is now easy to ob- 
tain for the bulk system by a Fourier transformation. Up 
to innocuous additive constants, the Hamiltonian of the 
vortex-free configuration now reads, 



where 



e q d\d q + i-£ (d\dl q + H.c, 



2J Z — 2J X cos q x — 2J y cos q y , 
2J X sin q x + 2J y sin q y . 



(19) 



(20) 



The fermionized Hamiltonian (|17p describes a p-wave 
type BCS pairing model with site-dependent chemical 
potential. After a Bogliubov transformation, this Hamil- 
tonian can be diagonalized and the quasiparticle excita- 
tion is 



(21) 



B. Spectrum 

The energy spectrum, in the low energy vortex-less sec- 
tor a r = 1 [and all sectors of {a r } related to it by the 
d = 1 operation of Eqs. (|13ll5[) ]. which we found to be 
given by Eq. ([20l l2Tj) obviously does not encode infor- 
mation about the topological nature of Kitaev's model. 
The equivalence of the lowest eigenvalues of the Hamil- 
tonians in both (i) a simple BCS type problem given by 
Eqs. (|2QI2ip and (ii) the Kitaev model of Eq.lU) vividly 
illustrates the maxim that the states themselves in a par- 
ticular (operator language) representation and not their 



q y 



J 7 



FIG. 4: A geometric interpretation of Eqs, (|20[l for the exis 



tence of gapless excitations. The two conditions e q 







imply that q XiV can be regarded as angles in the triangle 
formed by {J x , J y , J z } as shown. The law of cosines gives 
Eqs. (|22|l . This plot immediately leads to the triangle inequal- 
ity condition of Eq, (|27[) . When the triangle inequality is vi- 
olated (the gapped phase), q XtV are imaginary and lead to a 
finite correlation length as we will discuss later [Eq. (|67|) ]. 



energies which determine whether or not TQO exists. 
Q, Q In other words, the same Hamiltonian in different 
representations which are related to one another by uni- 
tary transformations and thus preserve the same set of 
eigenvalues (the energy spectrum) Q can describe both 
topologically quantum ordered systems (such as Kitaev's 
model) or systems with no topological order. Q Q In the 
current context, these mappings are the Jordan- Wigncr 
transformations that we apply. Q 

Let us now proceed to study when the spectrum of 
Eqs. pOl [2"Tj) describes a system with a spectral gap be- 
tween the ground and the next excited states and de- 
lineate these boundaries between the gapped and gapless 
phases. Towards this end- see Fig.Q- we study the nom- 
inal solution to E a = which is 



± arccos 



± arccos 




(22) 



This solution makes sense only when the conditions 



-2|J„||J*| < Jy 



-2\J X \\J Z \ < Jl + Jl 



Jl-J 2 x <2\J y \\J z 



Jy — 2 | J x | | J Z 



(23) 
(24) 



are satisfied simultaneously. These inequalities can be 
easily solved to obtain the condition for gapless excita- 
tions 



\Jx\ < \Jy\ + \Jz 

\jy\ < \ J X \ + \J Z 

\J Z \ < \J X \ + \J Z 



(25) 
(26) 
(27) 
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This result is in agreement with Kitaev's original solu- 
tion and Pachos' analysis. Q @ We can now view q xy 
as (up to an inversion q a — > — q a for one of the compo- 
nents of q) as the angles in the triangle formed by the 
couplings J Xl J y and J z - see Fig. (TJJ. Eq. ([22|) is the law 
of cosines in this triangle which automatically satisfies 
A g = e q = of Eq. ([20|) and consequently has E q = 0. It 
is also interesting to notice that the cyclic symmetry be- 
tween J x ,J y , J z is explicitly restored in the ground state 
solution although its explicitness is lost in the fermion- 
ized model. 

Next, we briefly remark on not only the low energy 
spectrum of pertinence to the zero temperature problem, 
but rather examine the entire spectrum of the theory. Af- 
ter tracing over the fermionic degrees of freedom present 
in Eq. (|17p . we obtain an effective two dimensional Ising 
type Hamiltonian (in {av}) with long range interactions. 
The full spectrum of this long range Ising type Hamilto- 
nian is identical to that of Kitaev's model. 



C. Real space form of the ground states 



The BCS type ground state corresponding to Eq. (|2T 



is 



\9) = n 



Uk 



Vk4dl 



10), 



(28) 



where 



E k 



1 + 2- 

E k 



(29) 



Let us now invert the transformations that we have per- 
formed until now (Jordan- Wigner and others) in order to 
express the Fermionic operators and the Fermionic vac- 
uum state in terms of the original spin degrees of free- 
dom. In what follows, we will, when needed, keep explicit 
track of the (w/b) sublattices of each of the sites of each 
vertical bond whose center is at f. Undoing all of the 
transformations that we employed thus far, we have 



4 = 



n 




[For a definition of e z see Fig.([TJ).] The product n r '<r 
corresponds to the product of all lattice sites f 1 which 
appear before r w on the Jordan- Wigner contour of Eq. © 
which traverses all sites of the two dimensional lattice. A 
similar definition applies to J] 



r'<r b 



it is the product over 
all lattice sites which appear before rf, in the Jordan- 
Wigncr product. 



In the spin basis, a fermionic vacuum corresponds to 




- 2 (l + h) 



(31) 



In Eq. (|3ip , the product over h is that over all elementary 
hexagons, \<j>) is an arbitrary reference state: e.g. in the 
a z basis, we may choose it to be the fully polarized state 
\4>) = | ft ... j*), the operator 



B r .„ = 



(32) 



extends over all sites lying between (and including) r w 
and rb as labeled by the one dimensional Jordan- Wigner 
contour that connects r w with ffe, J\f is a normalization 
factor, and the hexagonal operator Ih as given by Eq.©. 
Similarly, inverting the fermionization carried earlier, we 
find that 



n 

l<r w 



n <* 

r w <l'<r b 



(33) 



is the spin representation of the operator a r of Eq.([7|). 
Eq. (|33|) defines an operator a rw for any white lattice site 
r w . In the last product in Eq. (|3ip . we have a product of 
i(l + a rm ) over one lattice site r w in every row j. The 
symbol r* w denotes the first leftmost white lattice site 
r w in the jth row. In combination with Eq. (|12p [valid 
for any plaquette], the conditions a r » = 1 and 1^ = 1 
ensure that a r = 1 for all r. A derivation of Eqs. lpm i 
is given in the appendix. It is noteworthy that 



[4,B r J = [I h J h ,] = [B r „,S, 
= [Ih,a r * ] = [B T 



'J = [h,d 



0. 



(34) 



{Ih} and {B Tm } lead to disjoint Z2 algebras. Thus, when 
combined, Eql $M HD Hi W[ ESJ ED El E3D give us the 
explicit form of the ground state wavefunctions for the 
S = 1/2 system defined by Eq.|T|. 



D. Comparison to ground states of other 
topological ordered systems 

The spin basis form of the ground state of this topolog- 
ically ordered system can be compared to that of other 
systems exhibiting topological quantum order. E.g. both 
the Rokhsar-Kivelson [RK] [ll| and the Kitaev toric code 
models [H (as well as Wen's plaquette model [2l| which is 
equivalent to Kitaev's toric code model) @ have ground 
states of the form 



(35) 



with M a normalization constant and the sum over i per- 
formed over all states \xi) which belong to a given topo- 
logical sector which is evident by the real space represen- 
tation (e.g. an even or odd number of dimcrs/positivc 
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bonds in the Rokhsar-Kivelson/Kitaev toric code models 
respectively). In both the RK and Kitaev toric code mod- 
els, these ground states may be expressed as a product 
of projection operators (related to a algebra) acting 
on a given reference state. For instance, in the Kitaev's 
toric code model on the square lattice [4] 



with^ = n<oSandB p = n, 
all ground states may be expressed as 



(36) 



as stated earlier, 



Ifftoriccode) = ^\Yl\0- + A ')j 



v) 1$, (37) 



with \(j>) a reference state. Similarly, the general Eq. ((35|) 
can be written as a projection of a given reference state 
over a sector of Hilbert space that has a chosen set 
of topological numbers (those corresponding to a given 
topological sector). As seen from our solution of Eq. (f!?8")) . 
the ground state of Kitaev's model on the hexagonal lat- 
tice (EqfT]) [f| is more complicated. This is so as it in- 
volves the fermionic operators of Eq. (|30]) . The Fermi 
vacuum of Eq. (|3ip indeed already has a form similar to 
Eg. (|3T|) but, as seen from Eq. pS]) , it needs to be acted on 
by Fermionic string operators. This leads to a far more 
nontrivial state. Written longhand in the original spin 
basis, we have that the fermion pair creation operator 



4£* = zEE' 



,ik(r w2 —r w i) 



* [<j n <')-<,( n < 

\ r'<r w i r'<n,i 

(< 2 ( n <-)-< 2 ( n <*) 

\ r"<r W 2 r"<r b 2 



(38) 



needs to be augmented. Here, in the product signs we 
made explicit that along the Jordan- Wigner contour that 
we employ here r& = r w + L. When Eq. (|38p is combined 
with Eqs. lpM I3TI we have a form for the real space 
spin states of the Kitaev model. 



E. Time reversal symmetry breaking 

Our expressions show that the state \g) of the open 
boundary condition system is generally not time reversal 
invariant for any size system- whether it contains an odd 
or even number of spins. In the expressions above for 
the ground state this is seen by taking cr° — > — a a for all 
components a = x,y, z under time reversal. If any of the 
indices in Eq. (|38[) are such that [r w \ — r W 2) or L are odd 






FIG. 5: The geometry of three hexagons that span a total of 
thirteen vertices. By Kramers' theorem, ground states of this 
system that contains an odd number of S = 1/2 spins must 
break time reversal invariance. 



then Eq. (|38p is not time reversal invariant. Similarly, the 
operator £> rro of Eq. (|3"!2|) contains (L + 1) spin operators 
and is not time reversal invariant for even L. A reference 
state \<f>) of Eq. (|3~lj) on an even size lattice can be chosen 
to be time reversal invariant. For the odd size lattice, 
the considerations are more immediate. 

As a concrete example, let us consider three hexagons 
that span a total of thirteen vertices as shown in FigJS] 
This example serves to illustrate the presence of time 
reversal symmetry breaking even without doing any cal- 
culations. It is impossible to write down a ground state of 
Kitaev's model of thirteen S=l/2 spins of this thirteen 
site "disk like geometry" that is non-degenerate. The 
same considerations apply also to instances in which we 
have various topologies. For instance, we can put the fig- 
ure above of three hexagons with a total of thirteen sites 
on a three-fold torus- that of three handles (the genus 
number g = 3). Here, there will be a hole at the center 
of each of the hexagons. One of the sites is common to 
three ordinary tori and three other sites are common to 
two ordinary tori. There is periodicity around each of 
the hexagons (hence the designation of g = 3) in this 
case. The thirteen site system of Fig.© on the plane 
corresponds to a system with open boundary conditions 
while if it is a g = 3 torus is corresponds to a system 
that has non-trivial periodicity around 3 hexagonal rings 
(along with three other trivial periodicity directions of 
period one). In both cases, as the Hamiltonian H is time 
reversal invariant (because T~ 1 S a: iT = —S a -i with T the 
time reversal operator and S a; i the a = x,y, z compo- 
nent of the spin operator at spin i and the Hamiltonian 
is quadratic in the spins) and there is an odd number of 
spins (n = 13), by Kramers' theorem, the ground state- 
and in fact any energy eigenstate- must be, at least, two 
fold degenerate. General considerations for constructing 
a torus of genus q from individual ordinary (g = 1) tori 
are outlined in [1 61 ] - We emphasize that this degeneracy 
is mandated for an odd size system by Kramers' theorem 
for any time reversal invariant Hamiltonian. Time rever- 



7 



sal symmetry breaking similarly occurs in extensions of 
Kitaev's model that allow for odd cycles [j| [l5[ (e.g. a 
triangular decoration of the lattice). 

In any half integer spin system in which the number 
of particles is odd, the ground states must, by Kramers' 
theorem, exhibit time reversal symmetry breaking. This 
conclusion may be fortified for general gapped systems 
with short range Hamiltonians. In such systems, there is 
an exponential clustering of correlations. [TtJ This clus- 
tering, in turn, implies that the ground state can, up to 
exponentially small corrections, [181 ] be written as a ma- 
trix product state constructed of low lying states on finite 
size blocks. The block size may, in principle, be taken ar- 
bitrarily large so long as it is smaller than the system. 
Consider a general Hamiltonan H = y]^ fiR with R la- 
belling different blocks and /ir a Hamiltonian that has its 
support on R. Spins can be shared by different blocks R. 
By Hastings' theorem [l8| , the ground state of H is up 
to exponential corrections a matrix product state. Let 
us now examine what this implies when combined with 
time reversal. If all blocks R are chosen to have an odd 
number of sites in the decomposition of any short range 
5 = 1/2 Hamiltonian H, then all states of the different 
blocks R will be, at least, doubly degenerate. This is be- 
cause of the (at least) two fold degeneracy of each level 
in the odd size blocks R implied by Kramers' theorem. 
We see here that the matrix product construction then 
implies that the ground state of large systems (also if 
they are of even size) must also be, at least, two- fold de- 
generate up to corrections which are exponentially small 
in the size of the system. This degeneracy is dictated by 
the time reversal non invariance of the local energy cigcn- 
states in each of the odd sized blocks R. The considera- 
tions above are general and apply to both systems with 
or without topological order. Topological order pertains 
to systems such as Kitaev's which donot display global 
symmetry breaking associated with local order parame- 
ters. [l[ The existence of degeneracy in Kitaev's model 
can also be seen by the non-commutativity of existant 
symmetries: this non-commutativity lies at the heart at 
Kitaev's inception of this model. Wc will detail symme- 
try considerations for this model in a later section. 



F. Boundary terms 

For infinite open systems, the previous discussion will 
be sufficient for the analysis of bulk properties. How- 
ever, in a closed system, boundary terms will lead to a 
topological dependency of the ground states. In this sub- 
section, wc shall consider the toridal geometry, which can 
be viewed as periodic boundary conditions along both di- 
rections. For the Jordan- Wigcr transformation defined in 
Eq.©, the periodic boundary condition along the verti- 
cal direction will has no effect since the phase terms of 
two nearest neighboring terms along x or y bonds cancel 
out and there is no phase term for the z bond coupling. 
For the periodic boundary condition, the boundary terms 



read 

Hboundary = \jx&l,2j cr L h ,2j + ^y <7 l,2j + l IJ L h .2j + l $9) 

3 

Here, Lh is the size of the sytem along the horizontal 
direction. After the Jordan- Wigner transformation (J2J), 
the coupling strength of the boundary terms aquires a 
phase 

Rboundary ,2jA Lxt2 j 
j 

+J y <f>(2j + l)A h2j+1 A L ^ 2j+1 ] (40) 
with the phase term given by 

= II <r ( 41 ) 

l<i<L h 

This phase does not commute with the Z 2 bond operator 
a on the z-bonds attached to j-th horizontal line. The 
phase factor cf>(j) reverses all a of the z-bonds attached 
to the sites on the j-the horizontal line. The boundary 
terms thus lift the degeneracy characterized by differ- 
ent choices of a bonds. The gauge freedom discussed in 
section IIVI is now fixed by the boundary terms and the 
ground state degeneracy acquires a topological depen- 
dence. On the other hand, the boundary term commutes 
with the plaquette quantity defined by Kitaev. There- 
fore, the discussion about the correlator in later section 
is still valid when boundary term is included. The mo- 
mentum components (k Xl k v ) employed in the solution 
of the previous subsections above are related by a 45 de- 
gree rotation to the discrete values (kh, k z ) for the system 
with periodic boundary conditions along the vertical and 
horizontal directions discussed here. 

We conclude with brief remarks about the symmetries 
of Eq. p3)l in the case of periodic boundary conditions 
and on a candidate trial state. The spin system of Eq.fl]) 
on the torus has the cj>(j) operators of Eq. (fr51 l4"Tj) as 
exact symmetries. On the torus, the symmetry of Eq. (|15|) 
corresponds to the product of spin operators along an 
entire toric cycle. These operators satisfy the following 
identity 



<t>(j)<t>ti + 1) = II Ih 



(42) 



In Eq. (|4"2"|) , S is a ribbon of width one on the torus which 
contains all hexagons (plaquettes) lying between the two 
consecutive rows j and j + 1. As the ground states of 
the Kitaev model are vortex free (Ih = 1 for all h), the 
righthand side of Eq. (|4"2")) is one. Thus, for all rows j, we 
have that cj>(j) = <fi(j + 1). A variational state for the 
periodic system is given by 

IVw) = aa]T | fla o = a/-( + 4>U)]) \9) (43) 

a' j 

with M a normalization constant, and {a'} all configura- 
tions, which up to the symmetry operations of Eqs.(|13[ 
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ITS"]) , correspond to the uniform sector a r = 1 for all bonds 
r. Here, \g a >) is the ground state in a given a' sector of 
the open boundary system. Such a state is related to the 
general uniform a r = 1 ground state which we derived 
earlier (the state \g) of Eqs.jH [29] ffl El El])) by the 
application of the symmetry operators <f>(J) of Eqs. l|15[ 
HTJ . The variational state of Eq. (|4"3"|) is an eigenstate of 
all of the symmetries {<f)(j)} of Eq. lfTSI B4j) [with a uni- 
form eigenvalue which is equal to one] . 



VI. ANYONS IN THE GAPPED PHASE 

We will now derive and study anyons in the gapped 
phase of Kitaev's model by relying directly on the 
Jordan- Wigncr transformation. The existence of anyons 
in this phase has long been recognized [f| and, due to 
the prospect of fault tolerant quantum computing, is one 
of the main motivations for studying this system. The 
current appendix illustrates how anyons these may be 
directly studied and derived within our framework. 

Besides the Bogliubov quasiparticles, other 
excitations- the vortex excitation illustrated in Fig[3Jb)- 
are also manifest. It is more interesting to study these 
excitations within the gapped state. In this section, 
we shall follow Kitaev's original argument Q and now 
demonstrate, in our fermionic representation, the any- 
onic nature of these vortex excitations. Let us start 
in the limit where J x 
fcrmionized Hamiltonian reads 



J y = 0. In this limit, the 



H = J z ^a r (2did r -l) 



(44) 



The ground state is thus 2 N / 2 degenerate, where N is 
number of sites. This degeneracy can also be under- 
stood in the original spin language. In this limit, all 
vertical bonds are disconnected. Let us consider a single 
bond, (JzO-fuf). Without loss of generality, let us con- 
sider J z > 0. There are two degenerate ground states 
\]{) and | It)- This degeneracy is related to the local £T 2 
symmetry which present for J X J V = 0. Without losing 
generality, let us set J x = 0. In this case, we can define a 
local unitary transformation which is just erf erf on one of 
the x-bonds. The effect of this operator is to transform 
J z to —J z on the two vertical bonds connected to the 
x-bond. If J x J y ^ 0, such transformation is impossible 
and we expect the degeneracy to be lifted. The lowest or- 
der contribution from J x and J y term is thus expected to 
be in the second order of J x J y , i.e., J x Jy/J z . This cor- 
responds to a fourth-order perturbation. Indeed, Kitaev 
has shown this is true in the language of spin operators Q. 
The non-constant effective Hamiltonian up to 4-th order 
is thus proportional to the product of the four horizontal 
bonds of a given plaquette, (erf erf )(&2 a 3)( a 6 a 5)( a 5 a %)- 
In the language of our Majorana fermions A, this is 



H { e f f oc {A 1 A 2 )(A 2 A 3 )(A 6 A 5 ){A 5 A A ) 
= (iArAe^iAsAi). 



Here, we have dropped the sublattice subscripts b and w 
for the A fields. If we introduce a dual spin on vertical 
bonds, 



= 2c#cL - 1, 



(46) 



then the effective Hamiltonian H e ff will be nothing but 
an Ising coupling between two neighboring spins along 

The to- 



horizontal direction oc Mf ^+1/2/^+2 j+i/2' 

tal perturbativc Hamiltonian thus reads 



97 z Jx'Jy 

za ^ 1 ~ 16J3 



+ const. (47) 



on the lattice of vertical bonds. The coefficient —1/16 is 
worked out by KitaevQ. As one would expect, the J x 
and J y terms lift the degeneracy and the ground state is 
the one with on = 1 in our aforementioned gauge fixing 
of the extra freedom in section IIVI The interesting ex- 
citation state that corresponds to the anyon studied by 
KitaevQ is thus the ground state of the Hamiltonian of 
the configuration sketched in Fig[3jb). This is exactly 
the same vortex state in plaquette p defined by KitaevQ 
where one conserved quantity / is reversed in the plaque- 
tte p. In this state, to save (positive) J z contributions 
to the energy, the effective spins {fii} on the cut where 
cti = —1 are also flipped and the resultant configura- 
tion is a domain wall structure at p. The energy penalty 
associated with this domain wall isJ x Jy/8J z . For conve- 
nience let us define 



+ 1/2) = ^_i J+ i/2^+i J+ i/2- (48) 
In terms original spin operators, / is, similar to Eq.@, 
I = utulalalalal, (49) 

An illustration is provided in the lower corner plaquette 
of Fig® 

We shall now demonstrate that the anyonic nature of 
the vortex excitation. A vortex on the plaquette centered 
at (i,j + 1/2) is characterized by two kinks on the same 
bond, one in the channel of a and another in the channel 
of /i 2 . There are four equivalent ways to create a vortex 
on the plaquette centered at (i, j + 1/2): 



P*(i,j + l/2) = ]JiA w B w (i',j + l) (50) 

i' >i 

P%(i,j + 1/2) = Y[iA w B w (i',j + l) (51) 

i' <i 

P b R (i,j + 1/2) = Y[iA b B b (i',j) (52) 

i' >i 

P b L (i,j + 1/2) = l[iAiB h {i',j). (53) 

i' <i 

The operators (P b R ) are related to P^ (P b L ) by 
(45) the gauge transformations encapsulated by Uf (£/f) of 
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FIG. 6: Loop operator to move a vortex along a loop. 



J 5 , 



R,L 
b.w 



= 1. Thcrc- 



^b w can either create a vortex or annihilate an 



Ea.(fl"3]). We also notice that ( 
fore 

existing vortex at (i, j + 1/2). A naive way to move an 
anyon vertically from (i, j + l/2) to (i, j + 5/2) is to apply 
P b L (i, j + 5/2)P^(i,j + 1/2). It turns out that it is more 
convenient to multiply two extra phase terms, the first 
one contains the multiplication of z — bonds centered at 
+ 3/2) with i' < i 



J[{iA b B b )(iA w B w )} 



z — bonds 



and the second one contains the multiplication of I cen- 
tered at (i', j+3/2) withi' < i, l\ i , <i I(i',j+3/2). After 
including these two phases, the shifting operator T y from 
(i, j + 1/2) to (i, j + 5/2) takes a simple form in terms of 
spin operators 



T y = a x (i,j+l)o-y(i,j + 2). 



(54) 



We are now ready to construct the loop operator that 
moves vortex around and study the statistic properties 
of vortex excitations. Let us first create a vortex on the 
plaquette centered at (i, j + 1/2) with P^(i,j + 1/2), 
which applies S z on all site-3 (white site at top-right 
corner) of the plaquettes (i 1 , j + l/2) with i' > i as shown 
in Fig[5J To move this vortex along a closed loop c, we 
can construct a a loop operator Lc by multiplying T y 
along vertical lines and a z along horizontal parts of the 
loop, as illustrated in Fig[51 The loop operator can be 
expressed by 



Lc 



n /(*',/+ 1/2) [m 



(t',3')ec, 

j+j=odd 



'w u b\l 



(55) 



/ 



\hec 



An explanation of the notation is due. Without losing 
generality, let us assume the anyon is located on an even 



row, j = 2n. Lc can be separated into two parts. The 
first part is the product of L that covers the plaquettes 
on the odd rows enclosed by the loop C. The second 
part is the product of ct^ct^ on the vertical bonds whose 
black site is enclosed inside the loop C. In the projected 
subspace on which the effective Hamiltonian lives, the 
vertical bond part is the identity 1. Therefore, the loop 
operator is just the product of / over the plaquettes on 
the odd rows enclosed by C, Lc = ( — l)" c where nc is 
the number of vortices that are located on the odd rows 
and inside the loop C. In another words, if we bring a 
vortex on odd row around another vortex on an even row, 
an extra minus sign is generated, while no such minus 
sign if wc bring a vortex on an even (odd) row around 
another vortex also on an even (odd) row. Putting all 
of the pieces together, we fleshed out, by building on 
the fermionic representation, the anyonic nature of the 
vortex excitations. 



VII. STRING AND BRANE TYPE 
CORRELATORS DICTATED BY LOCAL 
SYMMETRY 

In this section, we shall apply symmetry considera- 
tions [H [IH and review condensed remarks from [l3[ 
concerning the symmetries of this model and general con- 
siderations regarding non-vanishing correlators in simi- 
lar matter coupled theories which are reviewed in the 
Appendix jB]). We will show directly from symmetry 
considerations that the Kitaev model on the hexagonal 
lattice [H does not exhibit any two point correlations of 
length larger than one and that any non-vanishing cor- 
relation function generally amounts to a string operator 
(cither closed or open). The fermionization that we em- 
ploy here [1] [lj| gives rise to precisely such string like 
correlators whose form is dictated by symmetry. In the 
matter coupled gauge analogy of Section jB]), all possi- 
ble correlators correspond to open meson or closed pho- 
ton lines. We will show that Kitaev's model supports 
finite valued brane type correlators. By brane correla- 
tors, wc allude to correlators which span portions of the 
two dimensional plane. This result- as it applies for only 
the two point correlations- was also recently derived by 
a Majorana fermion construction fl2| . In what follows, 
we show how gauge symmetry considerations effortlessly 
mandate this result. 



A. Two point correlators 

We claim that both at finite and at zero temperature 
the spin-spin correlation between any two spins which are 
separated by more than one lattice constant must vanish. 
To see this, we note that this system displays the local 
symmetries of Eq.®. We now fuse these symmetries 
together and see that we have a symmetry associated 



10 



with any closed contour C . by 

O c =l[<Tp pC . (56) 
pec 

In Eq. (|56p . the polarization index 7 p c is chosen to cor- 
respond to the direction (7 = x, y, z) of the single bond 
emanating from site p which does not lie in the contour 
C. Any given site p forms the endpoint of three bonds; 
two of these bonds lie in the contour C and only a single 
bond ending in p does not lie in C. 

We now briefly relate these symmetries to degenera- 
cies. If two curves C and C share an odd number of sites 
in common then Oc and Oc 1 anticommute. Thus if, e.g, 
we choose to simultaneously diagonalize H and Oq and 
determine a ground state |<?i) in that common eigenba- 
sis, then Oc \gi) is a new degenerate ground state. This 
degeneracy applies not only to the ground state sector 
but to all energy levels. 



1. T > 

By em ploy ing these symmetries along with Elitzur's 
theorem [lj| at temperatures T > 0, we find that the 
finite temperature correlator {(J^a b q ) vanishes unless (1) 
p and q are nearest neighbor sites and that (2) a = b 
corresponds to the direction (a = x, y, z) between the 
two lattice sites i and j in question. This is so as unless 
conditions (1) and (2) are satisfied, there is at least one 
loop C for which [0^0^] is not invariant under Oc of Eq. 
([56]) and for which 

6 c [o a p e\\6c = -[a a p a b q ]. (57) 

2. T=0 

The T = relation (even for different times [l2j) is sim- 
ilarly derived. All of the operators O c commute with one 
another. [This is seen as for the minimal loop C which 
contains a single hexagon, two neighboring hexagons 
leading to operators Oi and O2 which contain two com- 
mon sites which lead to two anti-commutation relations 
between the different Pauli matrices. [f|. The product 
of operators each acting of a single hexagon leads to the 
most general symmetry operator of Eq. (|56[) with an ar- 
bitrary path C] Consequently, we may simultaneously 
diagonalize H of Eq. (TTJ) with all of the operators O c . 
The dynamical T = two spin correlation function 

«K(o)o*(t)M 

= (VK(°) exp(iHt)a b q (0) ex P (-iHt)\ip) 

+^(^[H,[H,a b j m + ...) (58) 



with any ground state which is a simultaneous eigen- 
state of H and Oc- The eigenvalues of Oc can only be 
±1. We can now show that unless p and q are linked by 
a single step along a direction a each of the expectation 
values in Eq. (|58[) vanishes. For example, for the first 
term in the final expression, we have 

(vk^iv) = mocKa b (6 c m 

= (i>\(d c ay q dcM = ~{1>Ko*\1>), (59) 

unless p and q which are linked by a single step along 
a direction a. Similarly, the second term in Eq. (|58j) 
is seen to vanish. Here, the commutator [H, a b ] leads 
to sums of bilinears which involve the lattice site j and 
a nearest neighbor site k. When multiplied by a® this 
leads an expression which is cubic in the spin operators. 
If all of the three lattice sites are different, we can choose 
a path C such that this expression changes sign under 
a unitary transformation corresponding to Oc- Similar 
higher order terms in Eq. (|58[) unless p and q which are 
linked by a single step along a direction a. This result 
reaffirms (from a purely symmetry point of view) the 
more detailed derivation of 121. 



B. Higher order string and brane type correlators 

We now turn to higher order spin correlations. A mo- 
ment's reflection reveals that invariance under the gen- 
eral local symmetries of Eq. © allows only for the string 
or brane operators (either closed or open) to be finite: 
in other words, all spins must form continuous clusters 
along lines (strings) or lie in a fragment of the two di- 
mensional plane. 

String and brane type correlators may involve, in the 
thermodynamic limit, an infinite number of fields and al- 
low different topologies than that of closed loops alone. 
The fields may be ordered along an open string (or col- 
lection of such strings) or in higher dimensions may in- 
volve sophisticated combinations of fields at all lattice 
sites which we term as brane type correlators. 

The best known example of such a string correlator is 
that in the spin 5 = 1 AKLT chain [22j [23| in which 

fc-i . 
(Sf( II j" Si )Sk)= g- (60) 

Eq. (|6"0"]) holds for arbitrarily far separated sites i and 
k. As this correlator is, asymptotically far larger than 
the usual spin-spin correlator, it is often said to capture a 
hidden order. We may construct string operators which 
are more sophisticated variants of the operator appear- 
ing in Eq. (|60|) . Apart from spin chains, [22| [23( [2J] such 
correlators also appear, amongst others, in doped Hub- 
bard chains and related system s. l25ll [20 ] [27] ] cold atom 
chains [28|, in spin leg ladders |29l]. and may be related 
to non-local constructs in Quantum Hall systems [30j . 
In doped Hubbard chains, the string correlator decays 
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asymptotically with distance but with a power which is 
smaller than that of the usual spin-spin correlators. To 
date, nearly all appearances of string correlators are con- 
fined to one dimensional or quasi-one dimensional sys- 
tems. The only known exceptions which we are aware of 
are 12511 and possible links to the density matrix construct 
of [30( for Quantum Hall systems. In what follows, we 
show that the Kitacv model is a rigorous two dimensional 
example of a system in which various string correlators 
can either attain a finite value or decay asymptotically 
in a slow manner (explicitly so for gapless systems). 

Let us turn back to Kitaev's model of Eq.([T]). The only 
gauge invariant (symmetry of Eq.©) quantities which, 
by Elitzur's theorem [l!| can attain a finite expectation 
value are given by 



Q- 



P r> I 4 £ I 



Q- 



J* 



1 2 D 



-9- 



-G- 



-Q- 



-Q- 



-e- 



FIG. 7: An open loop correlator of the form given by Eg . (|61fl . 
Here, the correlator is equal ((erf erf )(<rf crfXcfcf)) = 
{erf a\ erfcrf). All correlators are of the form of continuous 



string correlators such as this or of their union. See Eq. (|61|l . 
As we review in the appendix, precisely such string correla- 
tors are the only non-vanishing correlation functions in matter 
coupled gauge theories. 



n 

\pgec 



Lp q 



(61) 



with C being a general contour which may be open or 
closed (or a union of such contours) and a pq denotes the 
direction of the physical bond between two nearest neigh- 
bor sites p and q. When C is a closed loop, the argument 
of the average in Eq. (|6ip is, up to a multiplicative phase 
factor, equal to the symmetry operator of Eq. ([56l) . 

An example of such a correlator is furnished in Fig. |(7J). 
In Fig.([7]l. we show the only non- vanishing correlator 
which has four spins with sites 1 and 4 as its endpoints. 
To see that this string correlator is invariant under all 
local symmetries, consider first sites 1 and 4. Consider- 
ing hexagons h = E, D (see figure), we see that only the 
fields erf and a% are invariant under the local symme- 
try of Eq.©. Next, we note that this particular choice 
of these fields at sites 1 and 4, the symmetry of Eq.© 
for both hexagons h = A and h = B, enables the in- 
troduction of a field at sites 2,3 of the form cr| 3 . Thus, 
{ (J \ (J \ <7 \ (J X) i s a correlator which is invariant under all 
local symmetries. As such, it is not prohibited from at- 
taining a non-zero expectation value at finite tempera- 
tures by Elitzur's theorem. Similarly, the product of two 
disjoint bonds ((erf c r f)( a x cr x)) ^ s invariant under all local 
symmetries. In a similar fashion, we can proceed to con- 
sider longer contours invariant under local symmetries- 
all of which must be of the form of Eg . (foTj) . Replicating 
the T = considerations of Eq. to multi-spin correla- 
tors, we see that string correlators of the form of Eq. (|6~lj) 
are similarly symmetry allowed within the ground state 
sector but others are not. By construction, the string cor- 
relator is invariant under all products of the local sym- 
metries \Jh\- the symmetries of Eq. (|56[) . The contour 
C may consist of disjoint open segments (e.g. disjoint 
bonds with each bond containing two sites). Of all of 
the correlators of Eq. fBTj) that we found to be allowed 
by Elitzur's theorem, only those with an even number 
of spins are time reversal invariant. String correlators 
of the form of Eq. (|6ip with an odd number of spins are 
not time reversal invariant and consequently must vanish 
when time reversal symmetry is unbroken. (2(j| 



For closed contours C, the quantity to be averaged 
in Eo.(|6ip becomes none other than the symmetry of 
Eq. (|56p . In Appendix (|B]), we review similar selection 
rules for lattice gauge theories. Though all non- vanishing 
correlators must be of the form of Eq. (|6"Tj) . systems with 
open boundary conditions, further allow for additional 
string correlators of this form with only the boundary 
site(s) on the contour C not adhering to this form. The 
appearance of only such non-vanishing open or closed 
continuous string correlators (involving all sites between 
sites p and q) is reminiscent of the gauge invariant cor- 
relators appearing in matter coupled gauge theories, e.g. 
[l3| . As stated earlier, all non- vanishing correlators must 
be of the form of Eq. (|6ip on a system with periodic 
boundary conditions or, in a system with open boundary 
conditions, further allow for additional string correlators 
of this form with only the boundary site(s) on the con- 
tour C not adhering to this form. The ideas underlying 
the constructs which we introduce next were reviewed in 



1. Correlators of maximal value 

We next construct several such string correlators by 
relying on the mapping to the Fermi problem. This map- 
ping transforms correlators which involve only several 
Fermi fields into those which involve an extensive number 
of spin fields and vice versa. 

In the Kitaev model, the Fermi vacuua states |0) can be 
chosen to simultaneously diagonalize all of the symmetry 
operators {Ih} and {B rm }. Here, we will have that 



; B r ,\o) 



l Cr r ral cr r rol +l Cr r rol +2 



< 3 -i<J0>| 



(62) 



(for an even number of sites). The string correlator of 
Eq. (|6"2"|) is equal to a constant value irrespective of its 
length (the number of spin operators that it contains). 
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Different Jordan- Wigner contours may also be chosen 
which lead to paths other than those in Eq. (f6"2")) which 
connect the two endpoints. The Fermi vacuua become 
ground states only in the limit of \J Z /Jx,y\ ^ 1 (and 
permutations thereof). Physically in this limit, the large 
value of \J Z \ forces consecutive chains to be fully corre- 
lated. The open string correlator of Eq. (|6"2"|) augments 
the closed loop correlators of Eq. ([SU)) for any closed con- 
tour C. Of course, combinations of various closed loops 
and open string operators are possible. 

We now discuss the situation for arbitrary J x ,y,z- 
Given a general ground state \g) of the form of Eq. (j2"8)) we 
can transform it into a Fermi vacuum state by a unitary 
transformation 

u = II ( Uk + Vk ( b l + 5fe ) - Uknk ) ' u \a) = |0>. (63) 



Here, n*, = d\.dk- For a general ground state |<?) the 
non-local brane expectation value 

|(.g|L/t< mi < ui+1 < iui+2 ...< b2 _ 1 < 2 t/| 5 )| = 1. (64) 

The operator of Eq. (f64|) has maximal correlations (the 
modulus of the expectation value of these operators can- 
not exceed 1). The maximal correlations exceed the stan- 
dard string correlators found in the AKLT chain. Writ- 
ten in terms of the original spin fields at each site, this 
operator generally (for finite J x ,y.z) spans all sites of the 
lattice. 



2. non-maximal correlators 

We can similarly construct other operators which 
would reduce to other open string operators if the system 
were one dimensional (if, e.g., J z =0). These form other 
generalizations of the familiar one dimensional string op- 
erators. In what follows, we discuss string operators 
which are not maximally (and indeed decay with increas- 
ing separation between the two endpoints of the string. 
To this end, let us start by writing 



4a 2 



±- 



1 



n 



<r lul <r'" <r wo 




n 



.(65) 



ir m , <r"<ri 



Here, the + sign is chosen if r w . 2 < along the contour 
on which the Jordan- Wigner transformation is performed 
and the — sign is chosen otherwise. Here, we have that 



1 



(4M - 2 / (27r)2 



d 2 k 



-ife-(ri-f 2 ) 



i-ii 

E k 



(66) 



with the definitions of Eq. (|21j) for and Ek- The form of 
this expression follows directly from the state of Eq. (f2"5)) . 
(rifc) = v\. Looking at Eq. (|66|) . we see that for gapless 
systems, a power law behavior may be sparked. In the 
presence of a gap, the branch points in the complex q 
plane of E q of Eq. (f2Tj) along a chosen direction may de- 
termine the asymptotic long distance correlation length 
along that direction. Rather explicitly, in the gapped 
phase, we have from Eq. (f22]) that the logarithm of the 
correlator of Eq. ([66|) scales asympotically, along the x,y 
directions as [-{\r 2 - a - ri; |/£ )] (where {r a } are r ■ e a 
with a = x,y, [see Fig.JT])]) with 



1 



cosh 



J: 



" z °y 



— = cosh 



1J X J Z 
J? + Jl - Jl 



^JyJ z 



(67) 



Along any given direction a, in the complex q a plane, 
there are two branch points along the imaginary q a axis. 
There is a simple signature of the T = critical phase 
that occurs when we cross from the gapped to the non- 
gapped phases when | J x ± J y \ = \ J Z \. Here, the branch 
points merge at q a = and there is a divergent correla- 
tion length £ a . Within the gapless phase, the two branch 
lie along the real q a axis. 

Fusing Eqs. (|65U66[) together, we have an expression for 
the expectation value of the string operator which lives 
on all sites linking r Wl up to (and including) . The left 
hand side of Eq. (|65[) is given by the Fourier transform in 
Eq. (|66p . Similar transformations may be applied to other 
fermionic correlators which when translated into the spin 
variables lead to other string correlators. 

The ground state form of Eq. (f2"8|) enables a direct com- 
putation of the expectation value of any general product 
of the form ( JT^ d/^ rj . d\ . } . This is so as the state of 
Eq. (|2"5|) is a direct product in k space. For each value of 
k, we have a binomial distribution for each occupancy, v\ 
for the occupancy of the pair (k, —k) and a probability of 
u'fc for it to be empty. Thus, any general multi-spin cor- 
relator can be trivially computed- we first Fourier trans- 
form it, then express the product in terms of the Fermi 
variables following a Jordan- Wigner transformation, and 
then compute the average for the states of Eq. l|2"8|) where 
each average for a given value of k becomes decoupled 
from all other values of k. In a related vein, Wick's the- 
orem for Eq. (|21[) ensures the decomposition of general 
correlator into a product of pair correlators for {k, —k) 
pairs. 



VIII. CONCLUSION 

In conclusion, we presented an exact solution of Ki- 
taev's toric code model which allows for various new re- 
sults: 

(1) We fermionize Kitaev's model to render it into a 
p-wave type BCS pairing problem. 
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(2) Wc derive the exact ground state wavefunction in 
terms of the spin degrees of freedom. With it, we exam- 
ine the zero temperature correlation functions and derive 
the exact correlation lengths within the gapped phase. 
We further show that the wavefunction in this case is 
more complicated than in the other prototpyical models 
of topological quantum order (e.g. the Rokhsar-Kivelson 
Quantum Dimcr Model and Kitaev's or Wen's square lat- 
tice models). 

(3) We prove that local symmetries only enable string 
(open or closed) and brane type correlators (and their 
unions) to be non-zero in this system. We show how 
string correlators may be directly evaluated in terms of 
the exact ground state wavefunction. Kitaev's model is 
one of the very few two dimensional systems with string 
or brane type correlators. Nearly all known examples of 
string correlators to date centered on one dimensional 
systems. 

(4) We identify the local symmetries of this system in 
terms of bond variables in the fermionic problem. 

(5) We illustrate, in terms of fermions, the anyonic 
character of the vortex excitations in the gapped phase 
by an explicit construction. 

Many possible extensions of our results follow from 
our fermionic wavefunctions. For instance, we can con- 
sider impurity bound states in systems engineered to have 
Hamiltonians which deviate slightly from Eq.|T]). For 
spatially non- uniform J x , J yi J z , so long as no new inter- 
actions are added, the toric symmetries remain in tact (as 
does the fermionization) . In the vortex free sector, the 
problem reduces to that of an impurity in a BCS type 
system which may lead to fermionic bound states. In 
subsequent work, we will further aim to detail the Ising 
type Hamiltonian that results for the vortex variables 
by tracing out over the fermionic degrees of freedom in 
Eq-CEZ). 
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Now, let us write the state |0) as a general superposition 
of states in the a z basis. 

|t7iCr 2 ...0-Ar). (A3) 

Inserting Eq. (|A2[) in Eq. (|A3|) we find that for all four 
possible spin orientations at r w and (up/down at each 
of these two sites) we have that the amplitudes satisfy 

Aai<n...crr m :.(Tr b :.(Tr N = — \ a r w +l ■ ■ -Or b -l ) 
A C r 1 cr2...cr rw -i(—a rw )a rw +i....<r rb -i(—cr Tb )a rb +i...try ■ (A-4) 

If and only if Eq. (|A4p is satisfied does the state of 
Eq.([A"3]) satisfy Eq.([A"2|). Eq. fAl)) is equivalent to the 
demand that 

«|0) = -« +1 ..<_ 1 )|0). (A5) 

In turn, Eq. (|A5[) is equivalent to the condition 

|0>=BrjO>, (A6) 

for all sites r w with the operator definition of Eq. (|3"2")) . 
Eq. (|A6[) along with the condition of no vortices as dic- 
tated by reflection positivity [f| (7/j = 1) or an easier 
immediate direct inspection of Eq. (fTTJ) leads to the gen- 
eral solution for the Fermi vacuum state in the spin basis 
fEq. (f3"Tj) ). It is noteworthy that in the particular limit 
J z > with J x = J y = 0,any of the 2 N / 2 states having 
a z ^a z h = 1 is a ground state. Here, N is the number of 
sites and (N/2) is the number of vertical [a z a z ) bonds. 
As seen from our fermionization (e.g. Eq. (|17p ). this cor- 
responds to the condition a r (2n r — 1) = —1 for each 
vertical bond r with n r = d\d r the fermionic occupancy. 
In the sector a r — 1 for all r, the remaining ground state 
is that of the fermionic vacuum derived above. As shown 
in Section pvp . the inversion of a r on all sites r of a given 
row leaves the spectrum invariant. If there are L horizon- 
tal rows then there are 2 L fermionic sectors that share 
the same spectrum. 

APPENDIX B: STRING CORRELATORS IN 
MATTER COUPLED GAUGE THEORIES 



APPENDIX A: THE FERMIONIC VACUUM IN 
THE SPIN REPRESENTATION 

In what follows, we briefly outline the derivation of 
Eq. (|3"Tj) . For the fermionic vacuum we have dk\0) =0 for 
all k or equivalently 



d r \0) = 



(Al) 



for all vertical bond centers r. Eq. (|Al[) implies that in 
the original spin basis, 



Kr'<r b 



l(T" 



n 



|0) = 0. (A2) 



Here, we briefly review, the well known local symme- 
tries of lattice gauge theories in order to clarify their simi- 
larity and the similarity of the string correlators that they 
mandate to the string correlators in Kitaev's model. In 
Section (jVTT|) . we invoke precisely this analogy. The crux 
of the selection rule on the allowed correlations (which 
forces all correlators be string type operators in both 
gauge theories and Kitaev's model) is Elitzur's theorem. 
Elitzur's theorem states that any quantity which does 
not transform as a singlet under local (gauge) symme- 
tries must have a vanishing expectation value. 

We now review Elitzur's theorem [l9| in its more 
prominent use- that of gauge theories. In theories of mat- 
ter at lattice sites (a p ) which are coupled to gauge fields 



14 



(U pq ) which reside on links of the lattice, [3l[ [1^ the ac- 
tion is a sum of (i) a plaquette product of the gauge fields 
and (ii) a minimal coupling between matter fields and the 
gauge fields. To illustrate, consider the simplest gauge 
theory [the Ising gauge theory] in which Ising gauge fields 
are coupled to Ising matter fields. Here, the action 

S = -if UpqUglUlnUnp - J ^ <J p U pq <J q , (Bl) 

with U pq = ±1 and a v = ±1. Many early results were 
found by [32| . This action is invariant under the local 
(gauge) transformations a p — > r/ p a p , U pq — ■> r] p U P qr] q with, 
at any site p, 7y p = ±1. Let us define gauge invariant 
link variables by z pq = a p U pq <j q . The action is a func- 
tional of {z P q}. Any correlator which involves a product 
of any number of z's is invariant under the local gauge 
transformations and consequently does need not vanish 
by Elitzur's theorem [lj|. In matter coupled gauge the- 



ories with J 7^ in Eq. (|Bl[) . any correlator of the form 

( n w (^2) 

pq ec 

for any contour C (either open or closed or a contour 
C which is the union of smaller open/closed contours) 
is gauge invariant and need not vanish by Elitzur's the- 
orem. For closed contours C, the average of Eq. (|B2p is 
the "Wilson loop" [3l| which makes an appearance also in 
"pure" ( J = 0) gauge theories. These string type correla- 
tors (as well as related topological percolation transitions 
and crossovers [33j) provide the only means for probing 
the behavior of this system. Precisely such string op- 
erators (both closed and open ended variants) are the 
sole non-zero correlators than are allowed by symmetry 
in Kitaev's model [see Section (|VIIBI) ]. 
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